Optimal feedback control for the rapid preparation of a single qubit 
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We consider the use of feedback control during a measurement to increase the rate at which a 
single qubit is purified, and more generally the rate at which near-pure states may be prepared. 
We derive the optimal bang-bang algorithm for rapid state preparation from an initially completely 
mixed state when the measurement basis is unrestricted, and evaluate its performance numerically. 
We also consider briefly the case in which the measurement basis is fixed with respect to the state 
to be prepared, and describe the qualitative structure of the optimal bang-bang algorithm. 
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I. INTRODUCTION 

In preparing a quantum system in a known pure state, 
one often starts with a mixed state, either because noise 
processes have rendered the initial state unknown, or be- 
cause the system is entangled with other quantum sys- 
tems. As a result, state-preparation requires the use of 
both measurement and unitary operations (Hamiltonian 
evolution). The measurement provides the required pu- 
rification, and the unitary evolution the ability to pick a 
particular final state. 

However, it turns out that for continuous measure- 
ments there exists a nontrivial interplay between the 
measurement process and unitary operations: the lat- 
ter, if applied during the measurement, change the rate 
at which the system is purified 0. Thus, a process of 
Hamiltonian feedback (that is, unitary operations per- 
formed conditionally upon the continuous output of the 
measurement process) plays the role not only of picking 
the appropriate final state, but also of determining the 
speed of projection, and thus ultimately of preparation. 
Note that while measurements are often treated as being 
instantaneous, all real measurements are continuous in 
that they take some time to act; measurements may be 
treated as instantaneous only if the measurement time- 
scale is much shorter than all other time-scales relevant 
to the given problem. 

The problem we consider here is the preparation of 
a single qubit in the fastest time. This is a particular 
instance of a problem in the domain of quantum feedback 
control 0, HH H H 0> and like a11 sucn problems it 
only makes sense when some or all of the resources at 
the disposal of the controller are finite. Thus we will be 
interested in the optimal speed of preparation given finite 
rates of measurement and/or Hamiltonian evolution. We 
will find that even this simple task of preparing a single 
qubit, employing the most straightforward measurement 
process, has a non-trivial structure. 

While the dynamics of closed quantum systems is lin- 
ear, that of observed quantum systems is often non- 
linear. The problem of feedback control of quantum sys- 
tems is therefore generally non-linear (even for very sim- 
ple systems), and our preparation problem is no excep- 
tion. Few analytic results exist for optimal control algo- 
rithms for non-linear systems. However, the tractability 



of the problem depends not only upon the dynamics of 
the system, but also the nature of the resource restric- 
tions and how they are applied. 

One method of enforcing resource restrictions is sim- 
ply to fix the maximal resources allowed, and to optimize 
the desired behavior of the algorithm under this restric- 
tion. In this case the functional to optimize is merely a 
function of the system dynamics under the action of the 
feedback. An alternative procedure is to make the re- 
sources flexible, and to minimize a 'cost' function which 
contains both a contribution from the amount of feedback 
required, and the resulting system dynamics. In this case 
'optimality' means optimality with respect to both the 
amount of resources used, and the resulting dynamics. 
The relative weighting between the cost of control and 
the cost of undesired dynamics is a free parameter. 

While the second of the above approaches is usually un- 
solvable for non- linear systems, a solution is often clear 
for the first if the constraint is merely on the strength of 
the applied feedback (the size of the applied force). That 
is, when one assumes that the rate at which this force 
can be changed is effectively unlimited on the relevant 
time scales of the system. The resulting algorithms are 
such that at any particular time the applied force is set 
at its extreme value in one direction or the other (for a 
one dimensional problem). The resulting feedback algo- 
rithms are referred to as 'bang-bang' control, a term that 
presumably originates from the sound that a mechanical 
controller makes when switching backwards and forwards 
between two extreme values. 

One can view our problem as consisting of a set of con- 
trol problems corresponding to increasingly realistic re- 
source constraints. We will find that when the strength 
of the measurement is fixed and the Hamiltonian is un- 
limited, the optimal algorithm may be found analytically, 
as well as its performance. When we next impose a limi- 
tation on the strength of the feedback Hamiltonian (but 
nothing else), there are broadly speaking two problems: 
one in which the measurement basis is free to vary and 
the second in which it is fixed. For the former the opti- 
mal bang-bang algorithm is easily obtained, but its per- 
formance must be calculated numerically. For the latter 
it is only possible to derive the overall structure of the al- 
gorithm. Both the quantitative details and performance 
must be calculated numerically. 
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In Section [H] we introduce continuous measurement 
and analyze rapid projection via feedback. Following 
reference Q we derive the optimal algorithm under the 
assumption of unlimited Hamiltonian resources, and its 
performance. The latter may be regarded as an an upper 
bound on all algorithms possible with a finite Hamilto- 
nian. In Section ITTTX1 we consider rapid state prepara- 
tion with finite Hamiltonian resources. In the first part 
we consider the case in which the measurement basis is 
under our control, derive the optimal algorithm and cal- 
culate numerically its perform as a function of the magni- 
tude of the Hamiltonian. This allows us to see how large 
a Hamiltonian is required to approach the upper bound 
in the previous section. In the second part we examine 
the case in which the measurement basis is fixed (and 
corresponds to the preparation basis), and describe the 
structure of the optimal bang-bang algorithm. 

Before we begin we note that recently there have ap- 
peared two related works by van Handel, Stockton and 
Mabuchi || |i| which consider the deterministic prepara- 
tion of one of the eigenstates of a continuously measured 
observable by using feedback during the measurement. 
There are also a number of articles which are concerned 
more generally with the feedback control of two-state 
quantum systems [3 13 13 13 El 

II. RAPID PURIFICATION VIA FEEDBACK 

All physical measurements are continuous - that is, ev- 
ery measurement extracts information at a finite rate. 
When this rate is much larger than all other relevant 
time-scales, then one can approximate the measurement 
by a von Neumann measurement. Otherwise one must 
describe the continuous extraction of information explic- 
itly. To describe a continuous measurement classically, 
one simply writes down the relationship of the measure- 
ment results, r(t) (which are a function of time), to the 
true value of the measured quantity, y(t), which is also 
in general a function of time. For virtually all applica- 
tions in which the stream of output results is a continuous 
function of time, this relationship is 

dr = y(t)dt + cdV, (1) 

where c is some constant, and dV is an increment of the 
Gaussian noise process referred to as the Wiener process. 
(For readers not familiar with the Wiener process, we 
note that easily accessible treatments may be found in 
references E3 and Gaussian noise is appropriate 

for virtually all continuous measurements in which r(t) a 
continuous function because of the central limit theorem. 
It turns out that the measurement record r may also be 
written in terms of the expectation value of y at each 
time. Specifically 

dr = (y(t))dt + cdW, (2) 

where dW is another zero-mean Wiener process (Gaus- 
sian white noise), uncorrelated with dV. 



One then uses Bayes' theorem 0, ^| to obtain an 
equation which gives the evolution of the observer's state- 
of- knowledge regarding y(t) (that is, her probability den- 
sity for y(t)) as the measurement progresses. This equa- 
tion is called the Kushner-Stratonovich equation E3 > and 
for a y which does not change with time, is given by 

dP(y, t) = 2V2k(y - (y))PdW, (3) 

where we have chosen k = l/(8c 2 ) to simplify the treat- 
ment in what follows. The quantity k characterizes 
the rate at which the measurement extracts informa- 
tion about y, and we will refer to it as the measurement 
strength. 

Quantum mechanically, one can describe a continuous 
measurement by using a sequence of "weak" measure- 
ments, where one takes the continuum limit of an infinite 
number of infinitely weak measurements pfjj . In doing so, 
one obtains the precise quantum equivalent of the classi- 
cal continuous measurement given by Eqs. J2Jl and |J3}. 
In the quantum case the observer's state-of-knowledge is 
given by the density matrix. If the measurement extracts 
information about an observable y, then the quantum 
equivalent of the Kushner-Stratonovich equation is 

dp = -k[y, [y, p]]dt + V2k{y P + py - 2{y)p)dW. (4) 

This equation is referred to as a Stochastic Master Equa- 
tion (SME). If y commutes with the initial density ma- 
trix, then this quantum equation reduces to the classical 
equation (Eq.©), because quantum measurement the- 
ory reduces to classical measurement theory (as it must) 
when all relevant quantities commute pi) . In this case 
P(y,t) is simply the diagonal of p(t). 

Let us consider now a continuous measurement of the 
z-component of spin of a spin-1/2 particle. The SME for 
this measurement is given by setting y = a z in Eq.QJ 
above. If we write the density matrix in terms of the 
Bloch vector, a = (a Xl a y , a z ), as 

p(i) = ~(/ + a-cx), (5) 

then the SME becomes 

da x = -(ikdt + a z V8kdW)a x , (6) 

da y = -(4kdt + a z V8kdW)a y , (7) 

da z = (1 - a 2 z )VskdW. (8) 

From these it is easy to show that the relation be- 
tween a x and a y is a constant of the motion, and thus 
the initial angle of the Bloch vector in the x-y plane, 
<p = aictan(a x / a y ) , remains unchanged throughout the 
measurement. We can thus reduce the equations of mo- 
tion to two variables. If we define A = ^a x + a y , being 

the length of the projection of the Bloch vector in the 
x-y plane, then the equations of motion become 

dA = -(4kdt + a z V8kdW)A, (9) 
da z = (1 - a 2 z )V8kdW. (10) 
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Due to symmetry the angle <f> plays no role in the dynam- 
ics. 

The purity of the density matrix, Tr[p 2 ], being the 
squared length of the Bloch vector, characterizes the ob- 
servers certainty regarding which pure state the system is 
in. If the observer has no information regarding the state 
of the system, this corresponds to a uniform distribution 
over all pure states (or, alternatively, equal probabilities 
for the two eigenstates of a z ). In this case the density 
matrix is described as being "completely mixed" , and is 
proportional to the identity matrix. In this case the pu- 
rity is minimal. If, on the other hand, the observer knows 
the system to be in a specific pure state, then the purity 
obtains its maximal value of unity. To measure the ob- 
servers uncertainty we will use p — 1 — Tr[/o], and refer to 
this as the impurity (the "bar" is intended to represent 
negation) . This quantity is often referred to as the linear 
entropy. In terms of the dynamical variables introduced 
above, the impurity is 



p=l-Tr[p 2 ] = l(l-A 2 



4). 



(11) 



If the state of the system is initially completely mixed 
(completely uncertain), then a 2 (0) = A(0) = 0. While in 
this case A remains zero throughout the measurement, 
a z (t) is stochastically driven to ±1. That is, the state 
is projected onto one of the eigenstates of tr z , which is 
the result of the measurement, and in this case each of 
these outcomes is selected with equal probability. In this 
case the measurement is also purely classical, since A 
is always zero, which means that the density matrix, p, 
commutes with o~ z at all times. 

The stochastic equations which describe the measure- 
ment process fEas. ljTU)) and ©) are non-linear. However, 
it is nevertheless possible to solve these equations by us- 
ing an equivalent formulation (referred to as a "linear 
quantum trajectory" ) in which they become linear 
(see also jUll). The solution is 



A(t) = csch(V8fcIF(i))A(0), (12) 

a z {t) = tavh(V8kW(t))a z (0), (13) 

where W(t) is a random variable, whose probability den- 
sity is 

p-ikt 

P{W, t) = -= cosh{V8kW)e~ w /(2t) (14) 
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^ e -(W+V8kt) 2 /(2t) +e -(W~Vskt) 2 /(2t)^ ( 15 ) 

Using this solution we can obtain an expression for the 
average impurity as the measurement proceeds. This is 



Pit) 



-4kt 



r-\-oc 



-x 2 /(2t) 



-dx. 



(16) 



cosh(v8fcx) 

While this integral must be solved numerically, certain 
features may nevertheless be extracted from this expres- 
sion for p. One finds that initially p decays as an expo- 
nential at rate 4fc, and an approximate analytic expres- 
sion may be derived in the long-time limit, to which we 
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FIG. 1: The maximum possible average speedup in the time 
required to achieve a final target purity from an initially com- 
pletely mixed state, for a single qubit, using Hamiltonian feed- 
back. This bound on speedup is obtainable in the limit of a 
large feedback Hamiltonian. 



will return later. In addition, general properties of mea- 
surement make it clear that p decays monotonically with 
time HIH. 

If this was a measurement on a classical bit, then the 
only way to increase the rate of reduction of the impurity 
would be to increase k, the strength of the measurement. 
However, quantum mechanically it is possible to increase 
this rate by using unitary operations during the measure- 
ment. To see this we must return to Eqs.© and I|1UI) . and 
calculate the rate of change of the square of the length 
of the Bloch vector, which is 



ri|a| ; 



riA 2 + dai 



(1 - a 2 )(A 2 - (1 + a 2 z ))8kdt 



.,[{1 - aj) - A 2 }V8kdW 



(17) 



It is clear from this equation that for a given value of |a| 2 
the average rate of increase is greatest when a z = (that 
is, when the Bloch vector lies in the x-y plane, perpen- 
dicular to the basis of the measurement). Furthermore, 
in this case the stochastic term vanishes, so that the in- 
crease in |a| 2 is deterministic. 

Since the rate of purification depends on the orienta- 
tion of the Bloch vector one can apply unitary operations 
to affect this rate (unitary operations can be used to ro- 
tate the Bloch vector, but cannot change the length). 
The unitary operator 

U(a) = exp[-i(a/2)(cos(0), - sin(</>), 0) • a], (18) 

will rotate the Bloch vector by an angle a towards or 
away from the x-y plane, while maintaining the angle 
(j). The Hamiltonian which generates this motion is thus 
H = hfj,(cos((j))a x — sm(4>)o~y), and the angle of rotation 
obtained is a = 2/j,t, where t is the duration over which 
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the Hamiltonian is in effect. The equations of motion 
describing such a rotation are simply 



dA = fia z dt, 
da z — —[iAdt. 



(19) 
(20) 



Now, we have seen from Eq. i|17|) that the rate of pu- 
rification is maximal when the Bloch vector lies in the 
x-y plane. However, from Eqs.® and it is clear 

that even if the Bloch vector lies in this plane at any 
given time, it will not remain there; the noise term in 
the equation of motion for a z will kick it out. (Note that 
if a z = 0, then the only stochastic term that remains 
is that in the equation for da z .) To maintain the most 
rapid average rate of purification we must therefore con- 
tinually rotate the Bloch vector back to the x-y plane as 
the measurement proceeds. If we allow the Hamiltonian 
to be arbitrarily large, then we can choose the rotation 
terms (EqsljSOl) to exactly cancel the noise term which 
kicks the Bloch vector out of the x-y plane. (The Hamil- 
tonian must be unbounded to achieve this, because dW 
scales as \/di) . Note that such a procedure is a feedback 
process, because the choice we make for the Hamiltonian 
at each point in time depends on the measurement result 
obtained at that time. In particular, to ensure that the 
Hamiltonian motion cancels the stochastic evolution due 
to the measurement, we must set 



/*(*) 



8fc dW 
~A~~dT' 



(21) 



There are a number of ways of calculating the evolution 
of the impurity under this feedback algorithm, but the 
most pedestrian turns out to be a little tricky; since the 
Hamiltonian term is now proportional to dW, we must 
consider the action of the Hamiltonian to second order. 
To perform the calculation we consider a single step of 
the measurement process, followed by a single step of the 
Hamiltonian feedback. At the start we set a z — 0, and at 
the end of the two steps a z is once again zero. The first 
step for A is given by Eq.@, with a z — 0, and adding 
the second step we have 



dA = -AkAdt + fia z dt - ^ 2 A(dt f, 



(22) 



where now a z — \/8kdW, being the a z which results from 
the first step, /i is given by Ea. (|21|l above, and we have 
included the action of the Hamiltonian to second order. 
Substituting in the relevant expressions, the result is 



dA = -4kdt{A - — ), 



(23) 



which is a purely deterministic differential equation. This 
might look a little odd - but if we write the equation for 
the impurity, p = 1/2 — 1/2A 2 , the behavior is immedi- 
ately clear: 



which is simply an exponential decay at rate 8fc. 

This result seems quite remarkable - not only does the 
Hamiltonian feedback increase the rate at which the sys- 
tem is projected, but greatly simplifies it. While the 
classical result (given by Ea. (|16fl % ) cannot be written in 
a fully analytic form, the quantum result with feedback 
is a simple exponential, and thus easily characterized by 
the rate 8fc. 

The feedback algorithm we have just presented is, in 
fact, optimal. The rate of increase of purity given by 
Ea. (|24|) is the best that can be achieved for any feed- 
back over any time period. While this is fairly clear from 
the above discussion, a systematic proof is given in ref- 
erence 0. 

So how much does the optimal feedback algorithm 
speed up the purification process? While we do not have 
an analytic form for the 'raw' purification as a function 
of time, we can obtain an expression valid in the limit 
as t — > oo. To do so we need merely note that the inte- 
gral in the expression for the raw impurity fEa. (|16ll ') is 
independent of t for large t. Thus for large t 



Pit) 



-4kt 



(25) 



for some constant C. If we set some very small impurity 
as the required target, then we can ask what is the ratio 
between the average classical time to reach this target, 
i c , and the optimal quantum time, t opt . When the target 
purity is sufficiently small (so that t c is sufficiently large), 
we have 



-4kt c 



or 



Ptargct 



*-opt 



In t c 

8kT c 



-8kt 



In 8irC 



8kt c 



(26) 



(27) 



dp = —Skpdt, 



(24) 



This tends to a value of 1/2 as t c — > oo, so that in the 
long-time limit the time taken to reach the target impu- 
rity is reduced by a factor of 2. A numerical solution of 
the integral expression for the raw p confirms that the 
speed-up factor increases with time monotonically, and 
thus the maximal speedup factor is 2. The speedup factor 
as a function of the target impurity is plotted in Figure^ 
While the optimal feedback algorithm presented above 
allows us to find the theoretical maximum extent to 
which unitary operations can increase the rate of en- 
tropy reduction rendered by a continuous measurement, 
the algorithm assumes unlimited Hamiltonian resources. 
We have not, therefore, determined how large a Hamil- 
tonian one would require in order to get close to the 
theoretical limit, and what speedup is achievable with 
a given strength of Hamiltonian. We have also been con- 
cerned merely with the purification of a quantum system, 
rather than the preparation of a particular pure quantum 
state. In the next section we will consider these ques- 
tions, and the implications of the above results for rapid 
state-preparation. 
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FIG. 2: The average speed-up in the time to obtain a purity 
of 0.95, from an initial purity of 5 x 10 -4 , for different values 
of the Hamiltonian feedback strength, A = Tr[iif 2 ]. The upper 
horizontal line is the maximal speed-up obtainable in the limit 
of large a large feedback, and the lower line is unity (no speed- 
up). Each plotted point was obtained by averaging over 4000 
realizations. Note that the plotted points do not fluctuate 
with respect to one another as much as the error bars would 
suggest - this is because the same random number sequence 
has been chosen in generating each point, with the result that 
the errors are highly correlated across the points. 



III. RAPID STATE-PREPARATION WITH 
FINITE HAMILTONIAN RESOURCES 



In the light of the above analysis, how can one prepare 
a given pure state in the fastest time? If the Hamiltonian 
resources at ones disposal are sufficiently large as to be 
effectively unlimited, then the answer is clear: one uses 
the Hamiltonian to keep the Bloch vector unbiased with 
respect to the measurement basis while the measurement 
proceeds, and then as soon as the desired level or purity 
has been achieved, one rotates the system to the desired 
state. However, what happens if the Hamiltonian re- 
sources are finite? To place a minimal restriction on the 
Hamiltonian we can put an upper bound on trace of its 
square: 



A. Flexible measurement 

If we wish to prepare a state \tp) from an initially com- 
pletely mixed state, and we are free to chose the basis in 
which our measurement is performed, then we can choose 
this basis to be unbiased with respect to the state \tp). 
Thus is equivalent to being able to choose our target state 
l^} as lying on the x-y plane, when we measure in the a z 
eigenbasis. 

The most rapid path to the state \ip) is the path which 
lies in the x-y plane, along the line on which the states 
have the same <fi as the target \ip)- (That is, the line 
which connects the center of the Bloch ball with the state 
This is clear because states in the x-y plane give the 
fastest rate of purification, and thus if (/) is chosen cor- 
rectly at the start (we will return to this point below), no 
further rotation is required at the end of the purification 
process. 

Since we have finite Hamiltonian resources, the feed- 
back algorithm will not be able to ensure that the state 
of the system remains exactly unbiased with respect to 
the measurement. However, since the average rate of 
purification decreases with increasing distance from the 
x-y plane, then, so long as 4> can be selected correctly at 
the start, the optimal algorithm is simply the one which 
keeps the state closest to this plane. Since we have not 
imposed restrictions on the rate at which our feedback 
can respond, the optimal algorithm is thus a bang-bang 
algorithm which rotates the state up or down towards the 
x-y plane with maximal speed, depending on whether the 
state is above or below the plane. The remaining issue 
is the initial selection of the angle 4>. At the start of the 
process, since we are assuming that the initial state is 
completely mixed, the initial action of the measurement 
is to point the state along the z-direction. Since this is 
the case, in choosing an initial Hamiltonian to rotate the 
state towards the x-y plane, we are also free to choose 
the initial (p. 

The optimal bang-bang algorithm for preparing a given 
nearly-pure state from an initially completely unknown 
state is thus clear. The dynamical equations, including 
this feedback, are now 



dA 

da z 

where 



-(4kdt + a z V8kdW)A + fi(t)a z dt, (29) 
(1 - a 2 z )V8kdW - (i(t)Adt. (30) 



Tr[iJ 2 ] < A 2 



(28) 



A , a z > 
-A , a z < 



(31) 



for some constant A. This restriction is unitarily invari- 
ant, so the controller is able to perform any rotation on 
the Bloch sphere, being limited only in the rate of this ro- 
tation. Under this resource limitation, one can identify 
two situations, one in which the basis of the measure- 
ment is free to move, and one it which it is fixed. We 
now consider each of these situations in turn. 



However, in contrast to the algorithm in the previous 
section, the resulting dynamics cannot be solved ex- 
actly; not even the technique of "linear quantum tra- 
jectories" [22, 13 i s sufficient. One must therefore solve 
these equations numerically. We do this and present the 
speed-up achieved for one value of the target purity as a 
function of feedback strength A in Figure [21 This shows 
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us that a feedback strength which is 100 times the mea- 
surement strength k is sufficient to get within 5% of the 
theoretical limit. 

We have examined the case in which the initial state 
of the qubit is completely mixed. The resulting opti- 
mal bang-bang algorithm is also applicable to the case 
in which the initial state is not completely mixed, but is 
a mixture of the identity and the state to be prepared - 
that is, lies along the line in the Bloch ball joining the 
target state to the center. However, this algorithm is no 
longer applicable when the initial state does not lie on 
this line. Further, it is not clear that the optimal algo- 
rithm in such a case can be obtained analytically. While 
we will not consider this problem in detail here, we de- 
scribe briefly the source of the complexity. When the 
initial state has the wrong </>, Hamiltonian resources are 
required to change the angle of the state in the x-y plane 
(that is, to change (j>). If the controller chooses to ro- 
tate the state in the x-y plane during the measurement, 
then this reduces the ability of the controller to keep the 
state in the x-y plane, and purification will be slower. 
However, if the controller waits until the desired purifi- 
cation has been achieved, and then rotates the state, this 
takes extra time at the end of the process. It is thus not 
clear how best the controller should apportion resources 
between purification and orientation during the measure- 
ment, and this remains an open problem for future work. 



B. Fixed measurement 

In the preceding section we considered rapid state 
preparation when the observer was able to freely choose 
the basis in which the measurement was made (i.e. the 
observable being measured). In this case we were able to 
obtain a simple optimal bang-bang algorithm for rapid 
state preparation when the initial state was completely 
mixed. In this section we consider the situation in which 
the observable is fixed with respect to the target state. In 
this case we will not be able to obtain the precise optimal 
bang-bang algorithm analytically, but we will be able to 
sketch the qualitative structure. 

If the measurement basis is fixed, then we cannot 
choose this basis to be unbiased with respect to the tar- 
get state. This means that the path of most rapid pu- 
rification will not include the target. Let us assume for 
concreteness that the measurement basis is the a z basis, 
and that the available Hamiltonian resources are such 
that the time required to rotate a state in the x-y plane 
to the target state is less than that required to obtain the 
target purity. If the controller decides to keep the state 
in the x-y plane until the required purity is reached, and 
then rotates the state to the target, the time required is 
the optimal time to purify, plus the final rotation time. 
It is clear that one can do better. This is because, the 
purification rate off the x-y plane, and thus closer to the 
target state, is not zero. Thus, if the controller starts to 
rotate the state towards the target shortly before obtain- 



ing the final purity, some time will be lost, but less than 
waiting until the end to perform the rotation. In fact, 
so long as the time taken to purify during the rotation 
is not greater than that required for the rotation, it is 
better to rotate before reaching the required purity. 

Conversely, a similar argument tells us that it is best 
not to rotate the state out of the x-y plane and towards 
the target if the average time to purification is greater 
than that required for the rotation. This is because, 
during any excess time required at the end for the re- 
maining purification, this purification will happen at a 
non-optimal rate. We can therefore conclude that there 
exists an optimal point of purity (distance from the cen- 
ter of the Bloch ball) at which it is no-longer desirable 
to keep the state on the x-y plane. Further, the same 
reasoning holds at any angle off the x-y plane towards 
the target state, but as the angle gets closer to the tar- 
get, the time required to rotate to the target is less. As 
a result, as the angle to the target decreases, the point 
of purity at which one should rotate towards the target 
is ever larger (ever closer to the edge of the Bloch ball), 
ultimately reaching the target purity at the target state. 

There exists therefore a critical surface, inside which 
one should rotate the state at the maximum available 
rate so as to keep it close to the the x-y plane. Outside 
this surface one should rotate the state as fast as possible 
towards the target state. That this rotation should be 
done as fast as possible follows from the fact that one 
wants to spend as little time as possible in the region 
off the x-y plane; one thus waits for as long as possible 
before rotating towards the target, and thus must rotate 
with maximum speed when the time comes. 

The qualitative form of the optimal bang-bang algo- 
rithm which results from the above reasoning is displayed 
in Figure 13 (for the case in which the measurement is in 
the a z basis, and the target state is |1)). The algorithm 
involves choosing the Hamiltonian so as to rotate the 
state towards the x-y plane (a z = 0), until a critical line 
is reached (the dotted lines in Figure [3] - there are two 
due to symmetry). These lines start at the x-y plane, 
and become closer to the edge of the Bloch ball as they 
move towards the target state, reaching the target purity 
at that point. Once the state of the system crosses these 
lines, the Hamiltonian is reversed so as to rotate towards 
the target. 

How might one go about calculating the exact form for 
the critical lines in Figure 13? The time taken to reach 
the target is ultimately a functional of the critical line. 
One could therefore simulate the feedback algorithm to 
obtain the average time-to-target, and seek to optimize 
this by varying the critical line. In implementing such 
a procedure one would ideally start with a guess for the 
critical curve based on approximate arguments. Even 
so, one would expect this method to require considerable 
numerical resources. 

The position of the critical line is determined by the 
average time required to reach the target purity from any 
given point. However, this average time is itself a func- 
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FIG. 3: A diagram giving the qualitative structure of the op- 
timal bang-bang feedback algorithm for the rapid preparation 
of a nearly-pure state close to |1), when the measurement is 
fixed in the z-basis {|0), |1)}. The grey portions denote ro- 
tation in the anti-clockwise direction at the maximum rate, 
and the white regions denote clockwise rotation. The pre- 
cise location of the dotted curves would have to be obtained 
numerically. 

tion of the feedback algorithm, and thus the positioning 
of the critical line - wherein lies the difficulty. It might 
be possible to obtain an approximate solution using a nu- 
merical algorithm which discretizes the angle to the tar- 
get, and considers small angles first. That is, calculates 
the remaining purification time at a small angle from the 



target, assuming that feedback at all larger angles in- 
volves rotating towards the target, and then once this is 
determined works backwards to larger angles. However, 
it is not clear that such an algorithm would reduce the 
numerical resource requirements. In any case, the task 
of obtaining the precise form of the optimal feedback al- 
gorithm remains an open problem for future work. 



IV. CONCLUSION 

We have seen that the rate at which a quantum state 
is purified by a continuous measurement is not isotropic 
on the Bloch sphere; it depends on the angle between 
the state and the measurement basis. This curious fact 
means that the task of rapid state preparation using feed- 
back control is non-trivial. It is possible to derive the 
optimal bang-bang algorithm for purification in a special 
case in which one is free to choose the measurement ba- 
sis, but it is not clear that a fully analytic solution will 
exist for the general case. Nevertheless, it is possible to 
describe the qualitative structure of the algorithm for the 
general case, and this algorithm does not involve taking 
the shortest path to the target state. 
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